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Abstract 

A  very  simple  derivation  of  the  Boltzraann  equation  is  given. 
It  is  based  on  the  solution  of  the  initial  value  problem  for  the 
"  hierarchy  equations"  of  statistical  mechanics,  i^iich  was  obtained 
In  research  report  No.  HT-3, 
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1,  Introduction 

In  an  earlier  paper^  -'  we  obtained  an  expansion  in  powers  of  the  density  of 
the  solution  of  the  initial  value  problem  for  the  "hierarchy  equations "  of 
statistical  mechanics.  The  purpose  of  the  present  paper  is  to  use  that  expansion 
to  obtain  a  very  simple  derivation  of  the  Boltzmann  equation. 

In  recent  years,  several  authors  have  given  derivations  of  the  Boltamann 
equation  based  on  the  hierarchy  equations.  J.  G.  Kirkwood*-  -'*'■•'  has  used  the 
hierarchy  equation  for  F,(the  one-particle  density  function),  and  by  means  of 
"phase  space  transformation  functions"  has  shown  that  the  Boltsmann  equation  is 
satisfied  by  a  function  obtained  from  F,  by  time-averaging.  N.  N.  Bogoliubov'-  -' 
has  obtained  the  Boltanann  equation  for  F,  by  assuming  that  for  a  >  1,  F  depends 
on  the  time  only  through  a  functional  dependence  on  F,.  M.  S.  Green ^■^■'  has 
obtained  the  Boltanann  equation  for  the  spatially  uniform  case  by  means  of 
"  Ursell  function"  expansions  of  solutions  of  the  hierarchy  equations. 

For  the  present  derivation  we  are  indebted  to  all  of  these,  but  our  task  is 
greatly  simplified  by  the  fact  that  we  have  at  our  disposal  an  explicit  solution 
of  the  initial  value  problem  for  the  hierarchy  equations.  The  leading  term  of  the 
expansion  of  that  solution  consists  of  an  integral  which  already  closely  resembles 
the  Boltmann  collision  integral.  Following  the  method  of  Green,  we  simplify  this 
integral  by  making  use  of  the  fact  that  in  all  but  a  small  portion  of  the  initial 
configuration  space  of  two  particles,  two-body  interactions  can  be  described  in 
tenns  of  complete  collisions.  Under  the  assumptions  usually  made  in  deriving  the 
Boltzmann  equation  we  find,  in  agreement  with  Kirkwood,  that  that  equation  is 
satisfied  by  a  time-averaged  density  function.  It  appears  likely  that  by  using 
further  terms  of  our  expansion,  which  involve  interactions  of  more  than  two  particles, 
the  method  used  here  can  be  extended  to  obtain  generalizations  of  the  Boltzmann 
aquation  to  higher  densities* 
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2.     Derivation  of  the  Boltzmann  Equation. 

The  statistical  mechanical  behavior  of  a  classical  system  of  identical 
monatomic  particles  can  be  described  by  the  infinite  set  of  equations 


<^>     ^-[V^]*i 
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d^3^lJ  3  -  1,   2, 


Here  F  (t,x.,  •••»3Cg)  is  the  s-partide  density  function,  t  denotes  the  time, 
X.  -  (q.*p J  -  (qyqj>qj,p.,p.>pi),   the  q.  and  p     are  the  cartesian  coordinates 

and  conjugate  momenta  of  the  J       partd.cle,  0  is  the  inter-particle  potential 
function,  V  Is  the  specific  rolunie,  and  H     is  the  hamlltonlan  for  a  system  of  s 
particles  t 


(2) 


y"T(p^)+Ug;  Ug  -V"  0(|qi-qj|)   J  T(p^) 

1=1  i<r<.i<8 


i<r<j<8 

The  mass  of  a  particle  Is  denoted  by  m.     In  (l)  Pbisson  bracket  notation  is  used 
and  the  integration  extends  over  the  entii«  x    , -space. 

In  [l]  we  obtained  an  esqianslon  of  the  solution  of  the  initial  value  problem 
for  (1)  in  the  fom 
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Here  S^"''  is  the  solution  operator  of  an  s-particle  systemi  i.e.  if  a  system  of 
t 

8  particles  (with  hamlltonlan  H  )  is  represented  at  time  t  -  0  by  the  state 
jx,,  ...,  X  I,  at  time  t  it  will  be  represented  by  the  state |x-,  .«.,  x  i-  si' 


jx.,  ...,  x\.     If  g  is  a  function  of  (z,  x^, 


.(s). 


X  ^  )  we  define  S^^g  by 
s+m  z 


These  equations  are  called  the  "  hierarchy  equations "  or  sometimes  the 
"B.B.O.K.T.  equations"  (for  Bogollubov,  Bona,  H.  S.  Green,  Kixkwood  and  Tvon). 


the  equation 
.(s) 


(u)  s;«^g(.,x^, ...,  x^^)  -  g(z,s;='^  |x^, ...,  x^  X3,i,  ...SX^J 


An  alternative  form  of  the  expansion  (3)  is  given  by 
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where  the  operator  T^°  is  defined  by 


.(s) 


s  -  1,  2,  ...} 


(1), 


,(1), 


(6)  Ti'^g(z,x^,  ...,x^j  -  g(z,s;«^  [x,,  ...,X3) ,  s;^\^i,  ...,s;^^x3^) . 


Formally,  (5)  can  be  obtained  from  (3)  by  first  interchanging  integration  and 

stmmation  over  the  index  J  and  then  introducing  the  transformation  x._^.  ^^  ■ 

S  i  xl+3+1,  ....  X  ^.  -  S^i^x'  .  of  the  integration  variables.  The  first  step  leads 

••w   J  3+K     "ti   S+K 

to  divergent  integrals.  However,  the  same  transformation  can  be  performed 
rigorously  at  a  conv^enient  point  in  the  derivation  of  (3), 

The  leading  terms  of  (5)  for  s  ■  1  are  given  by 

(7)  F^(t,x^)  -  T^i\(o,x^)  *  I  I  JTifF2(o,x^,X2)-T^i^F2(o,x^,X2}ldX2  *   O(^)  . 

Since  initial  data  may  be  specified  at  an  arbitrary  time,  (?)  may  be  rewritten 
as 

(8)  T^i%*V,x^)  -  "riWit^x^)  *^  j   feS^*'^'^2^  -"^^t^  F2(t,x^,X2)  dXg  ♦  O(^)   . 

■It 

Thus  S^''  acts  on  the  first  s  of  the  variables  x.  appearing  in  g. 

The  transfomation  may  be  carried  out  f.n  the  ecjuation  preceeding  equation  (38)  of    [l] , 


Since  (8)  is  an  identity  in  x,  we  may  replace  x,  by  s\-  x, .  If  we  then  introduce 
the  transformation  of  the  integration  variable  Xj  »  s!^  ■'x'  we  obtain 

(9)  i[^(t*L,S<,l'x,).Fi(t.:=^)]  .  i  j  [s(2Vj(t^(l>Xi,s(,^)xj)-Fj(t..^,Xj)]dqjdPj«)(^ 

We  now  introduce  a  time-average  of  F,  defined  by 

(10)  F^(t,x^)  -ij  F^(t+s,  S^-'-^)  ds  . 

We  shall  show  that  F^  satisfies  the  Baltzmann  equation  under  the  following 
assumptions t 

(I)  0(r)  ~  0        for  r  >  r.    . 

(n)  FgCt.x^jXg)  «   F^(t,Xj^)F^(t,X2)     for  hg-q^^l  >  r^,  for  some  r^  >  0   . 

(Ill)  T  »  1  . 

(17)  Fj^(t,q^  ♦  /^q^^,  p^^)  «   F^(t,q^,Pj^)     where   \/\fl^\  »  t^,t^     . 

(II)  is  a  generalized  form  of  the  molecular  chaos  assumption.     This  form  of 
the  assumption  and  the  general  technique  which  will  be  used  to  obtain  the 
collision  integral  were  suggested  by  a  paper  of  M.  S.  Green '■•^-'.     A  detailed 
discussion  of  the  asstimption,  which  is  part  of  Green's  "  product  condition"  is 
contained  in  [3].     In  (17)  we  assume  that  F.  does  not  change  appreciably  under 
translations  AQi  where   l/^q^  I  is  large  compared  to  both  r,(the  range  of 
intennolecular  forces)     and  t^     (the  correlation  distance)   .       F^  may, 
however,   change  considerably  over  distances  of  the  order  of  a  mean  free  path. 
Conditions  (I)  and  (ni)  are  self-cxplanatoiy.   (I),   (II),  (III),  and  (17) 
are  the  assumptions  usually  made  in  deriving  the  Boltzmarn  equation.     The  earlier 


The  operation  that  leads  from  F,  to  T.   is  sometimes  called   ■  coarse-graining" 

or  "  smoothing"  .     The  latter  term  would  appear  to  be  preferable.     See   [2],  part  I. 


deidvations  reqxiired  that  (II)  hold  for  all  values  of  |q«-q^|. 
Let 


a^        ^<'l 


Then,  since  S^  ^x^  -  S^     (q^^jP^)  •  (^i*^  Pj^jP^)  we  have 


(i2)Vi(W-sV**"=iS> 


2-0      f 


i  J^F^(t.s,S^^>x^)ds  -i[F^(t*t;,S;;^l5x^).F^(t,X3L)]   . 


We  recognize  (12)  as  the  leftsdde  of  (9)  • 

In  order  to  evaluate  the  right  side  of  (9)  we  choose  fixed  values  of  q^,  p., 
and  p«  and  examine  the  integration  with  respect  to  q^.  The  region  of  integration 
may  be  conveniently  analyzed  as  in  Fig.  1. 


Figure  I 
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A,  B,  and  C  together  form  a  cylindrical  region  with  spherical  end  surfaces. 
Outside  of  this  region  S_^  {s!j?"^x^,S^-'"'x2l  -  {3Cj^»3C2  f  *"^  *^  integrand  in  (9) 
vanishes.  For  points  in  C,  the  operator  s}^      mapo  q^   into  q^  (the  position  of 
S^^'qg  relative  to  fixed  q^^)  and  S_^^  maps  Qg  into  q^  .  For  all  points  in  C 

having  the  same  orthogonal  projection  (given  by  the  vector  b)  onto  the  plane 

t      I  t      I 

P,  the  momenta  p.  and  p  are  the  same,  i.e.  p^  and  p-  are  fVinctions  of  p.,  p. 

and  the  "  Impact  vector"  b.  From  (9)  and  (12)  we  have 

(13)  D^F3^(t,qpP^) -i  J  dpg  j  dA  F^(t,qpP^)F^(t,qpPp-F^(t,q^,Pj^)Fi(t,q;^,P2)l|-~— I 

*ii.o(^). 


The  integral  in  (13)  represents  the  contribution  to  the  integral  in  (9)  from  the 
region  C  except  for  an  error  due  to  the  fact  that  the  ends  of  the  cylinder  C  are 
not  planes  but  sphertcal  segments.  This  error,  plus  the  contributions  frcsn  the 
regions  A  and  B,  determine  the  tenn;^  I  .  In  (13)  we  have  clearly  made  use  of 
assumptions  (n)  and  (17).  The  volvune  element  dqg  has  been  replaced  by -Ipg-p^^ldA. 
where  dA  is  an  element  of  area  on  the  plane  P.  Now  from  (9) 

(lU)  F^(t+8,S^^^x^)  -  F^(t,x^)  ♦  0(|). 

By  Inserting  (Ub)  in  (10)  we  obtain 

(15)  F^_(t,X3^)  -  F3^(t,x^)  +  0(|)  . 

This  result  enables  us  to  replace  F^^  by  f^   in  (13),  the  only  change  being  in  the 
term  of  order  -^. 


-  7  - 

Let  u  be  the  average  relative  velocity  of  two  particles.  If  we  now  choose 
t   such  that  tu  »   r.(j»l,2),  the  region  C  will  be  much  larger  than  the  regions 

A  and  B  and  the  term  —-I  may  be  neglected.  By  (III)  we  may  omit  the  term  of 

order  -k  and  the  resulting  equation  is 

V 


(16)  D^F^(t,qj^,p^)-  I   dpg  j  dAk(t,q3L,pi)F^(t,q^,P2)-yj^(t,q,p^)T'^(t,qj^,P2)l 


^1. 


This  equation  is  the  Boltzmann  equation. 
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